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CHAPTER  1.  BASIC  CONCEPTS  AND  METHODOLOGY 


Section  1 . 1  Introduction 

The  concept  of  availability  represents  the  desire  that  a  piece  of 
equipment  be  operational  when  it  is  called  upon  to  function.  Availability 
extends  the  basic  notion  of  reliability  which  in  its  pure  meaning  only 
deals  with  the  continuing  proper  operation  of  a  unit.  Availability  allows 
equipment  to  undergo  repair  or  replacement  either  due  to  breakdown  or 
as  standard  maintainance  procedures.  The  concept  of  availability  can 
be  applied  to  a  small  unit,  complex  systems,  or  generalized  to  personnel. 
It  therefore  forms  a  major  part  in  logistic  analyses  over  a  broad  area  of 
applications.  A  number  of  mathematical  models  have  been  developed  for 
the  analysis  and  design  of  availability  requirements.  An  example  of 
such  mathematical  approaches  to  some  aspects  of  availability  are  given 
in  (1)  and  (2)  where  a  definition  is  developed  for  component  availability 
which,  unlike  component  reliability,  must  be  placed  within  the  context 
of  the  system  in  which  a  component  is  operating. 

The  field  of  contract  structuring  is  also  concerned  with  availability. 
Most  contracts  for  the  procurement  of  complex  systems  include  specifi¬ 
cation  of  the  availability  of  an  operational  system  for  a  period  following 
acceptance  of  the  equipment  itself.  One  approach  to  this  aspect  of 
contract  structuring  has  been  through  the  Reliability  Improvement 
Warranty  (RIW).  By  requiring  the  contractor  to  assume  some  form  of 
responsiblity  for  equipment  failure,  the  RIW  gives  a  cost  measure  to 
reliability.  The  usual  measurement  of  reliability  in  terms  of  probability 
concepts  like  mean  time  between  failures  is  difficult  to  deal  with  in 
contract  negotiations,  acceptance  tests,  and  disputes  over  peformance. 
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By  introducing  a  cost  measure  for  reliability,  the  concept  becomes  more 
clearly  recognized  as  a  real  feature  of  contract  responsibility  and  its 
dimensions  can  be  spelled  out  in  specific  terms.  Of  course  what  is 
really  desired  is  a  warranty  on  availability,  which  includes  both  reli¬ 
ability  concepts  and  aspects  of  the  logistic  support  required  for  equip¬ 
ment  supplied  under  contract.  In  (3)  and  (4)  the  ideas  of  availability 
warranties  and  incentive  type  contract  structures  have  been  given 
mathematical  representation.  That  material  illustrates  the  mathematical 
model  approach  to  the  analysis  of  availability  warranties. 

Most  of  the  mathematical  models  and  analyses  of  availability  and 
availability  (reliability)  warranties  have  considered  individual  units  of 
equipment  or  collections  of  such  units  in  which  the  availability  of  each 
unit  is  desirable  and  important.  One  may  refer  to  such  considerations 
as  "unit  availability."  The  major  mathematical  tools  for  such  studies  are 
probability,  statistics,  and  queuing  theory.  An  illustration  of  this  kind 
of  analysis  is  given  in  (5).  The  advantage  of  probabilistic  models  is 
that  a  deterministic  analysis  of  unit  failure  is  in  most  cases  too  complex 
and  the  methodology  of  stochastic  processes  is  well  developed  and  yields 
useful  results  under  many  conditions.  However,  stochastic  models  have 
two  important  disadvantages:  they  require  data  for  their  useful  imple¬ 
mentation  that  are  often  very  difficult  or  indeed  impossible  to  obtain, 
and  they  are  hard  to  interpret  to  people  untrained  in  the  details  of 
probabilistic  reasoning,  for  example  as  part  of  a  contract  negotiation. 
Thus  it  would  be  desirable  to  have  alternative  methodologies  for  avail¬ 
ability  models  if  they  could  be  developed  so  as  to  give  meaningful  and 
useful  analysis.  This  does  not  seem  possible  in  most  cases.  However, 
for  the  concept  of  bulk  availability,  such  an  alternative  methodology  does 


seem  possible.  This  is  the  primary  goal  of  the  work  described  in  the 
report:  to  define  a  non-probabilistic  model  for  bulk  availability.  It  is 
felt  that  development  of  such  a  model  will  show  the  feasibility  of  an 
alternative  methodology  for  at  least  some  types  of  availability  analyses. 
It  also  provides  a  model  that  is  useful  directly  for  the  study  of  bulk 
availability.  In  particular,  it  gives  a  basis  for  the  formulation  of 
availability  warranties  in  contracts  that  deal  with  bulk  availability. 

Bulk  availability  differs  from  unit  availabiltiy  in  that  the  individual 
units  are  not  the  primary  consideration  but  rather  the  number  of  units 
that  are  available  at  a  given  time.  A  bulk  availability  model  applies  to 
situations  in  which  there  are  a  number  of  similar  units  making  up  the 
total  system  which  is  defined  as  the  collection  of  units.  Thus  the 
system  itself  is  unstructured  in  the  model  formulation.  Availability  of 
such  a  bulk  system  is  defined  in  terms  of  the  units  that  are  available. 
Such  models  apply  to  a  group  of  workers  that  are  essentially  inter¬ 
changeable  for  functional  purposes  such  as  a  platoon  of  army  personnel. 
The  platoon  is  operational  if  some  specified  percent  of  its  composite 
personnel  is  operational.  In  the  same  way  the  model  applies  to  fleets  of 
trucks  or  aircraft,  highly  redundant  groups  of  equipment,  or  gracefully 
degrading  complex  equipment. 

Because  of  the  definition  of  bulk  availability  as  a  fraction  of  the 
total  collection  of  units  being  operational  a  deterministic  model  suggests 
itself  as  an  alternative  to  the  usual  stochastic  models  of  availability.  In 
bulk  availability  the  individual  units  are  not  important  but  only  the 
number  that  are  in  a  particular  state  of  operation  at  a  particular  time. 
The  model  presented  in  this  report  employs  transfer  rates  between 
states  such  that  the  actual  number  of  units  changing  state  depend  on 


these  rates  and  the  number  of  units  in  a  state.  Other  model  definitions 
are  also  possible  and  may  be  more  appropriate  in  some  situations. 
However,  the  present  model,  as  defined  in  the  next  section,  applies  to 
situations  in  which  the  number  of  transfers  from  one  state  to  another 
depend  on  the  number  of  units  in  a  state  and  the  transfer  rate.  For 
example  in  a  platoon  that  is  exposed  to  a  sickness,  the  number  of 
personnel  who  become  ill  in  each  time  period  depends  on  the  rate  of 
contagion  and  on  the  number  still  unaffected  who  are  therefore  able  to 
become  ill. 

Section  6  employs  the  bulk  availability  model  to  describe  a  method¬ 
ology  for  the  analysis  of  bulk  availability  warranties  and  related  logistic 
considerations.  This  material  is  based  on  (6)  while  this  report  gives 
some  additional  background  and  extensions  of  the  basic  ideas  given  in 
(6). 

Section  1.2  Definitions  and  Terminology 

This  study  considers  a  model  for  bulk  availability  in  which  there 
are  a  large  number  of  similar  units  that  fail  at  a  constant  rate  Upon 
failure,  a  unit  enters  a  waiting  line  and  from  there  enters  repair  service. 
The  service  activity  can  accomodate  up  to  r  units  and  completes  service 
at  a  constant  rate  p.  It  is  assumed  that  the  number  of  units  departing 
from  a  state  is  equal  to  the  appropriate  rate  parameter  times  the  number 
of  units  required  or  available  depending  on  which  state  is  under  con¬ 
sideration.  The  model  is  constructed  as  a  system  of  ordinary  differential 
equations  for  the  three  quantities: 

rn  (t)  =  the  number  of  units  that  are  active  (available), 

a 

mw(t)  =  the  number  of  units  that  are  waiting  for  service,  and 
m  (t)  =  the  number  of  units  in  service. 
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All  three  quantities  are  functions  of  time  as  the  independent  variable. 

In  this  model  system,  availability  can  be  defined  as  a  lower  bound 
on  rri  (t). 

a 

It  is  assumed  that  the  system  is  conservative  in  that  ma(t)  +  mw(t) 
+  ms(t)  is  a  constant. 

The  differential  equation  system  is  intrinsically  non-linear  in  that 
part  of  the  forcing  function  depends  on  the  values  of  the  unknown 
functions  for  its  formulation.  Thus  the  model  is  formulated  in  two 
regions  of  the  (m  ,m  )  plane.  Figure  1  shows  the  relations  between 

o  W 

the  three  states  of  the  system:  available,  waiting,  and  in  service. 


no  size  limit 


no  size  limit 


at  most  r 

Service  /u.  ms(t) 


Figure  1 

States  of  the  Bulk  Availability  Model 


The  system  of  differential  equations  governing  the  flow  of  units 
through  the  states  of  the  bulk  availability  model  follows  directly  from 
the  conservation  and  flow  rate  assumptions  stated  above.  This  system 
has  the  form: 


dma(t) 


dmw(t) 


dms(t) 


Mms(t)  -  Ama(t) 


'\ma(t)  -  K(t) 


=  K(t)  -  um„(t) 


In  this  system,  the  quantity  K(t)  depends  on  the  relative  values 


of  mg  and  mw  leading  to  two  specific  forms  for  K(t)  as  given  below. 
Because  of  the  conservation  assumption,  the  initial  values  sx  =  ms(tQ), 
vvi  =  anc*  ai  =  m*(0  satisfy  the  condition  s*  +  wx  +  a!  =  m  (t) 

+  mw(t)  +  ma(t)  for  all  values  of  t.  Since  the  differential  equation 
system  must  be  followed  across  different  forms  the  "initial"  values 
a ! , w L ,  and  st  are  incorporated  as  parameters  into  the  solutions. 

The  value  of  K(t)  depends  on  what  is  required  by  the  service 
state  and  what  is  available  from  the  waiting  state.  This  leads  to  the 
following  division  of  an  (mg,  mw)  plane  into  two  major  regions. 

Region  1  K  =  m  (t) ,  defined  by  m  £  r  and  m  <  r  -  m  +  pm  . 

VY  O  W  b  b 

Region  2  K  =  r  -  m  (t)  +  pm  (t),  defined  by  m  S  r  and 

b  b  S 

mw  =  r  -  ms  +  Mms- 
These  regions  are  shown  in  Figure  2. 


mW 


Solutions  are  obtained  for  Regions  1  and  2  in  Sections  3  and  4 
respectively.  Some  numerical  examples  are  given  in  Section  5  to  illus¬ 
trate  bulk  availability  analysis  using  the  model  described  here. 


One  procedure  for  using  the  model  would  be  to  simulate  the  system 
using  the  differential  equations  to  provide  the  model  structure.  This  is 
equivalent  to  direct  numerical  solution  of  the  equations.  However,  the 
system  of  differential  equations  can  be  solved  directly  in  terms  of 
explicit  functions  of  time.  This  is  the  approach  taken  here  and  given 
in  the  following  sections.  The  resulting  functions  of  time  are  compli¬ 
cated  and  require  computer  evaluation  as  described  in  Section  5.  Even 
so,  direct  explicit  solutions  are  felt  to  be  superior  to  numerical  solu¬ 
tions  of  the  differential  equation  system.  In  particular  such  explicit 
solutions  can  be  studied  directly  for  limiting  values  and  special  cases  of 
the  model  parameters. 

Section  1.3  Scope  of  Model  Methodology 

The  deterministic  mathematical  model  for  bulk  availability  given  in 
this  paper  relates  the  availability  of  a  system  to  various  system  para¬ 
meters.  Unlike  the  usual  stochastic  model  it  does  not  depend  upon  an 
underlying  probability  framework  and  its  results  are  specific  quantities 
expressed  as  functions  of  time  rather  than  expected  values  or  other 
partial  descriptions  of  random  events.  In  this  sense,  the  present  bulk 
availability  model  is  related  to  stochastic  models  in  the  way  thermody¬ 
namics  is  related  to  statistical  mechanics  as  a  description  of  certain 
physical  processes.  The  goal  of  thermodynamics  is  to  relate  various 
measurable  quantities,  e.g.  temperature,  to  physical  parameters  without 
the  necessity  of  a  detailed  analysis  of  the  partical  dynamics  that 
produce  such  relations.  This  is  also  true  of  the  bulk  availability  model 
in  which  the  transition  rates  and  initial  distribution  of  units  among 
states  are  system  parameters.  The  model  relates  these  parameters  to 


the  descriptive  functions  m  ,m  ,  and  m  . 
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The  value  of  the  model  is  in  its  direct  relation  between  a  few 
system  parameters  and  clearly  identified  descriptors  of  the  system, 
particularly  rn  (t)  which  measures  system  availability.  This  direct 

a 

relation  is  useful  in  contract  negotiations,  system  design,  logistic  sup¬ 
port  design  and  implementation,  and  the  conducting  of  acceptance  test 
procedures.  Stochastic  models  are  much  harder  to  interpret  for  use  in 
all  these  aspects  of  contracting  activity.  However,  this  valuable 
feature  of  the  bulk  availability  model  is  obtained  at  the  expense  of 
rather  limited  applicability.  For  the  model  to  be  meaningful,  the 
assumptions  required  for  the  model  must  hold.  This  restricts  model 
applicability  to  the  bulk  case  in  which  the  system  description  can 
indeed  be  given  by  a  division  of  units  into  the  three  states  of  the 
model.  Thus  the  model  takes  no  account  of  any  interaction  between 
units  or  differences  in  utility  between  units.  The  conservation 
assumption  also  limits  the  model.  There  are  only  the  three  specified 
states  so  that  any  case  in  which  a  unit  fails  and  cannot  be  repaired  or 
replaced  is  not  included  in  the  present  form  of  the  model.  Of  course 
such  additional  states  can  be  included  in  the  model.  The  desirability  of 
doing  so  to  extend  model  validity  must  be  balanced  against  the  desire 
for  relatively  simple  system  descriptors  which  motivated  this  kind  of 
model  from  the  beginning. 

The  most  serious  limitation  of  the  model  is  the  assumption  that  the 
number  of  units  that  transfer  state  can  be  expressed  as  a  product  of  a 
transfer  rate  and  the  amount  available  for  transfer.  In  any  particular 
case,  this  may  or  may  not  be  a  valid  assumption.  Its  validity  can  be 
established  by  theoretical  arguments  or  assumptions  or  by  experimental 
data.  If  it  cannot  be  validated  in  some  way,  it  is  likely  that  the  par- 
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ticular  bulk  availability  model  presented  here  does  not  apply  and  should 
not  be  used.  Other  assumptions  can  .  t  made  in  a  deterministic  model 
context,  but  the  ones  used  here  seem  to  give  the  most  direct  example 
of  the  kind  of  model  of  interest  in  this  study. 

CHAPTER  2.  REGION  1  SOLUTIONS 

Section  2 . 1  General  Form  of  Region  1  Solutions 

As  specified  in  Section  1.2,  the  Region  1  form  of  the  bulk  avail¬ 
ability  model  is  given  by  the  value  K(t)  =  mw(t).  This  results  in  the 
following  system  of  differential  equations: 

ma(t)  =  Mms(t)  -  Ama(t), 

mw(t)  =  Ama(t)  -  mw(t), 

ms(t)  =  mw(t)  -  nms(t), 

where  the  dot  notation  is  used  to  indicated  differentiation  with  respect 
to  time.  The  initial  conditions  for  the  system  are  values  ma(tQ),  mw(t0), 
and  ms(tQ)  specified  at  a  time  t  which  represents  the  time  at  which 
this  system  starts  to  govern  the  behavior  of  the  model.  Over  time,  the 
model  may  change  from  a  Region  1  to  a  Region  2  system  or  conversely. 
Therefore,  the  initial  conditions  for  a  system,  including  the  initial  time 
may  be  obtained  as  specified  data  or  as  the  final  solutions  of  a  previous 
solution  form.  This  is  described  more  fully  in  Section  4  on  numerical 
solution  procedures. 

The  method  of  Laplace  transforms  is  used  to  obtain  the  Region  1 
solutions.  If  L(g(t)]  =  G(s)  is  the  Laplace  transform  of  a  function  g(t), 
define:  L[ma(t)]  =  Fa(s),  L[mw(t)]  =  Fw(s),  and  L[mg(t)]  =  Fg(s). 


Also  define  the  initial  values  as:  ma(tQ)  =  ai<  row(tQ)  =  Wl'  anc* 
m  (t  )  =  Si-  The  transformed  equations  for  the  Region  1  system  have 

b  U 

the  following  standard  form: 

(s+A)Fa  -  pFs  =  ax 

-\Fa  +  (s+l)Fw  =  Wl 

-Fw  +  (s+M)Fs  =  Sl 


This  transformed  system  can  be  solved  in  the  form: 


=  ^F 


s+A  s  s+A 


w 


Am  r,  .  3l  wi 

(s+l)(s+A)  Fs  +  (s+l)(s+A)  ^  +  1+1 


where 


s^fs+lXs+A)  w,(s+A)  a,  A 

P(s)  P(s)  P(s)' 


and  P(s)  =  sQ(s),  where  Q(s)  =  s2  +  (1+A+m)s  +  A  +  m  +  M*-- 


The  form  of  the  solution  depends  on  the  roots  of  the  quadratic  function 
Q(s).  The  descriminant  for  this  quadratic  is: 

D  =  (1-A)2  -  2m(1+A)  +  m2 
or  alternatively: 

D  =  (1-A-m)2  -  4mA. 

The  curve  D  =  0  is  a  rotated  parabola  symmetric  about  the  line 


A  =  m  in  the  (m,A)  plane  as  shown  in  Figure  3. 
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Figure  3 

Areas  of  (p,A)  Plane  Defining  D  Values 


A  more  detailed  plot  of  the  curve  D  =  0  is  shown  in  Figure  4. 
The  region  with  D  <  0  is  inside  the  parabola  as  shown.  Any  combi¬ 
nation  of  p  and  A,  within  this  region,  gives  a  Case  1  solution.  Values 
outside  the  parabola,  where  D  >  0,  correspond  to  Case  3  solutions.  In 
the  model  as  presently  formulated  p  and  A  assume  fixed  values  in  an 
analysis.  However,  in  generalizing  the  model  or  in  considering  the 
kind  of  results  an  analysis  would  yield  it  is  interesting  to  consider 


linear  variation  of  the  form  p  =  KA.  Straight  lines  such  as  these  inter¬ 
sect  the  parabola  in  two  points  given  by  the  A  values  A^  =  ■ 

When  K  =  1,  there  is  only  one  solution,  K  >  1  gives  solutions  on  the 


lower  branch  and  K  <  1  gives  solutions  on  the  upper  branch  of  the 


parabola  respectively.  A  typical  intersection  is  shown  in  Figure  4  for 
the  value  K  =  2  for  which  A^  =  3±2j2  giving  values  of  (.344,  .172)  and 
(11.656,  5.828)  as  the  two  points  of  intersection  shown.  The  behavior 


of  the  solutions  changes  as  p  varies  along  a  curve  p  =  KA.  For  small 


values  of  p,  a  Case  3  type  solution  occurs.  As  p  increased  the  first 
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H- 


Figure  4 

Regions  of  the  (n,\)  Plane 


intersection  point  is  reached,  where  D  =  0,  giving  a  Case  2  solution. 
Further  increase  in  p  moves  the  solution  through  a  range  where  D  <  0 
and  a  Case  1  type  solution  governs  the  model.  For  K  =  1,  this  situ¬ 
ation  continues.  For  all  other  values  of  K,  increase  in  p  takes  the 
solution  to  the  second  intersection,  a  Case  2  form,  then  on  to  Case  3 
form  solutions  again.  The  regions  specified  in  Figure  4  can  be  helpful 
in  determining  the  general  form  of  solutions  to  be  expected  when  using 
various  (p,\)  combinations  and  to  see  how  the  solutions  will  change  as 
the  (p,A.)  values  change. 

The  roots  of  Q(s)  are  given  by:  -a±J  ^  where  a  =  (l+\+p)/2. 
Let  p2  =  D/4.  Then  the  values  of  D  as  a  function  of  p  and  \  yield 
three  distinct  forms  of  solution  as  indicated  by  the  three  regions  of  the 
^p^£)  plane  shown  in  Figure  3.  These  cases  are  defined  by:  D  <  0 
gives  Case  1,  D  =  0  gives  Case  2,  and  D  >  0  gives  Case  3. 

In  Case  1  the  quadratic  Q(s)  =  (s+a)2  +  p2;  in  Case  2, 
Q(s)  =  (s+a)2;  and  in  Case  3,  Q(s)  =  (s+a)2  -  p2.  The  form  of  the 
roots  of  Q(s)  in  terms  of  p  can  be  treated  in  various  ways  and  care 
must  be  exercised  to  get  the  correct  form  of  the  solution  in  each  of  the 
cases  and  also  in  numerical  calculations  as  discussed  in  Section  4 


Section  2.2  Case  1  Solutions 

In  this  case,  D  <  0  and  Q(s)  has  complex  roots.  The  quantity  p2 
is  negative.  It  is  useful  to  represent  Q(s)  as  a  quadratic  of  the  form 
Q(s)  =  (s+a)2  -  p2.  Then  Fg  can  be  expressed  as  follows: 


F  =  SjAfs)  +  wxB(s)  +  ajC(s),  where 


A(s)  =  s2  1  (1+A.)s  +  A 
s[(s+a)2  -  p2] 


At 


BjS  +  C  i 
(s+a)2  -  p2 


s 
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Direct  calculation  yields: 


Ai 


X 

X+|J+X|J  ' 


Bi 


U+Xm  r  _  |j+Xm+X2p 

X+p+Xp  '  1  X+p+Xp 


In  some  formulas,  it  is  convenient  to  express  X+p+Xjj  as  a2  - 
which  for  Case  1  gives  a2  -  p2  since  p2  =  D/4.  However,  it  must  be 
observed  that  in  this  case,  a2  -  p2  is  in  fact  the  sum  of  two  positive 
quantities  in  order  to  get  the  correct  inverse  Laplace  transform  expres¬ 
sion.  Proceeding,  we  find 


B(s)  = 


s  +  X 

s[(s+a)2  -  p2] 


A2  B2S  +  C2 

—  +  -  ' 

s  (s+a)2  -  p1 


where 

A2  =  Ai  =  X+p+Xp  '  ®2  =  X+p+Xp  ' 


c  = 

2  X+p+Xp 


C(s)  =  - - - - 

s[(s+a)2  -  p2] 


A3  B3S  +  C3 

-  +  -  , 

s  (s+a)2  -  p2 


where  A3  =  A1(  B3  =  B2, 


_  -X(a+X+p) 
X+|J+X|J 


The  inverse  transforms  are: 


-1  -at  (-Cl  "  Bl°^  -at  . 

L  i[A(s)]=  Ai  +  B te  cospt  +  - p - e  smpt. 


-1  -at  (C2  ”  B2a)  -at  . 

L  i[B(s)]=  Aj  +  B2e  a  cospt  +  - p - e  UL  sinpt, 


-1  -at  ('Cs  ”  B2°^  -at  . 

L  i(C(s)]=  At  +  B2e  cospt  +  - p - e  UL  sinpt, 


m_(t)  =  s1L‘1[A(s)]  +  w1L"1[B(s>]  +  aiL  1[C(s)]. 


Thus: 


^ID 
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Inversion  of  this  transform  yields  m,  ,(t).  The  Case  1  solutions,  ob- 

w 

tained  by  means  of  Laplace  transforms  as  shown  above,  are: 

1  _/yt 

m  (t)  =  AA  +  -  [(jsx(l+\)  -  Ka^w^Je  cospt 

s  a2-p2 

+  - — -  [(jSjCl+A2-^-^)  +  Wi(2p-A.2+A+A(j) 

2p(a2-p2) 

-  \a1(l+\+|j)  Je"0't  sinpt 

rn  (t)  =  pA  +  — - —  [-(Si+w^M  +  a1\(l+(j)]e"C!t  cospt 
a  a2-p2 

+  - - -  [s1jj(A+|j+2Aij-1)  -  Wj|j(l+\+|j) 

2p(a2-p2) 

+  a1A(l+p2-A-A(j)]e~0lt  sinpt 

rn  7(t)  =  pAA  +  — - —  [w1(A+(j)  -  Ap(s1+a1)]e  cospt 
w  a2-p2 

+  - 1 -  [w1(A2+p2-A-p)  +  a^CZA+p+Ap-p2) 

2P(«2-P2) 

«/yt 

-  sxAp(l+p+A)]e  sinpt 
In  these  equations,  A  =  Ax. 

As  a  check  on  these  rather  involved  results,  it  is  true  that 
ma(t)  +  mw(t)  +  mg(t)  =  aL  +  w:  +  sx  as  required  by  the  conservation 
of  units  assumption  in  the  model  formulation.  This  direct  though  some¬ 
what  tedious  calculation  is  omitted  from  this  report.  Some  additional 
checks  on  the  calculations  will  be  discussed  at  the  end  of  the  chapter. 
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Section  2.3  Case  2  Solutions 

In  Case  2,  D  =  0  so  that  Q(s)  has  two  real  roots  both  equal  to  -a. 
Thus  Q(s)  =  (s+a)2  so  that  the  transform  of  rng(t)  can  be  written  as: 

Sis2  +  [s1(l+\)  +  wjs  +  (Si+Wx+ax)A 

F  =  -  , 

s(s+a)2 


=  A(s)  +  (Sx+Wx+a^ABCs), 


SxS  +  Sx(l+A)  +  wx  , 

where  A(s)  =  -  ,  B(s)  =  - 

(s+a)2  s(s+a)2 


Partial  fraction  expansion  yields: 

w(  +  Sx(l+\-a)  Sx 

A(s)  =  - ; -  + 

(s+a)2  s  ( 


s+a  ' 


B(s)  =  - — - - - - = -  ,  so  that: 

a2s  a(s+a)2  a2(s+a) 


(s1+w1+a1)A  .  (sx+wx+ax)^  , 

Fs  =  - I -  I  +  ls> - j -  1 

oh  oh 


(Sx+Wx+ax)A  , 

+  [wj  +  Sx(l+A-a) - - - ]  - 

(s+a)5 


The  transform  of  mw(t)  is  given  by  Fw 
above  expression  for  F  ,  we  can  write 


(s+m)F’s  "  si' 


From  the 


-  +  R2  sL  +  Rs 


(s+a)5 


where  R  x  = 


(Sx+Wx+ai)A 


r2  =  Si  -  Rx, 


R3  =  Wx  +  Sx(l+A-a)  -  aRx- 
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Then  F, ,  can  be  written 
w 

Fw  *  R.  ^  *  R3-^  -  S 

(s+a)2 

Partial  fraction  expansion  and  simplification  yields: 


1 


Fw  ‘  i  +  [R2(M-«)+R3]^  + 

The  transform  of  ma(t)  is  given  by: 


ai 

Fa  =  +  MRxA(s)  +  pR2B(s)  +  pR3C(s), 

A(s)  =  sTs^T)  '  B(s)  =  (s+XXs+a)  '  and 
Partial  fraction  expansion  yields: 

A(s)  =  7a  "  XT^X)  '  B(S)  =  (a-X^s+X)  '  T< 

and  C(s)  =  - - - - - - -  . 

(a-X)2(s+X)  (a-X)(s+a)2  (a-X)2(s+a) 

These  quantities  give  F_  in  the  following  form: 

a 


mR 

Fa  X 


m^i  i  pRi  mR2  mR3  i 

X  s  +  X  +  a-X  +  (a-X)2  ^  s+X 


♦  [- 


mR2  mR3 


a-X 


(a-X)2 


1  JL  ^  1 

S+a  "  a’K  (s+a)5 


'  •  •  1 

The  coefficient  of  is  equal  to  zero  as  can 
direct  algebraic  calculation. 


where 

C(s)  =  - ± - 

(s+X)(s+a)' 

1 

-X)(s+a)  ' 


be  established  by 
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Inversion  of  the  transforms  developed  above  yields  the  Case  2 
solutions.  To  simplify  the  notation  in  these  solutions,  let  B  = 
(Si+Wj+aiVa2,  and  C  =  Wi+s^l+A-a)  -  AorB.  Then  the  Case  2  solutions 
are: 


m  (t)  =  AB  +  [sx-AB]e'at  +  Cte"at 


m  (t)  =  pB  +  — *= -  [-s1-w1+A(2a-A)B]e 

3  (A-a)2 


-at  .  Cte-cft 
of-p 


mir(t)  =  pAB  +  [wx  -  ApB]e  +  (p-a)Cte  at 
w 

By  direct  calculation,  it  can  be  shown  that  the  conservation 
assumption  is  satisfied  for  these  solution  functions,  i.e.,  that 
m  (t)+m  (t)+m  (t)  =  a^Wj+Si  for  all  values  of  t. 

a  W  O 

Section  2.4  Case  3  Solution 

In  this  case,  D  >  0  and  Q(s)  has  two  distinct  real  roots: 
9X  =  -a  +  02  =  -a  -  It  is  convenient  to  sometimes  express  ^ 

as  p  in  the  solution  forms  for  Case  3.  In  this  notation,  the  roots  of 
Q(s)  are  -a+p,  and  -a-p. 

The  transform  of  ms(t)  in  Case  3  has  the  form: 


s,(s+l)(s+A) 


Wi (s+A) 


s(s-e1)(s-e2)  s(s-e1)(s-e2) 

s(s-e1)(s-e2) 

A(s)  +  (a1+w1+s1)AB(s), 

sxs  +  Sx(l+A)  +  Wj 

1 

Expansion  gives: 


A('s-)  S“6 !  +  s-e2  * 


where 


S101  +  Sj(l+\)  +  wt 


S  i02  +  Sx(l+X)  +  Wx 


a2  b2  p, 

B(s)  =  ~s  +  ~s^q[  +  -i^e7  - 


where  A2  = 


0i02  ' 


B2  = 


0i(0i-e2) 


and  C  = 


02(0i-02> 


These  expressions  result  in  the  value: 


Hxi+  H2  -4—  +  H3-4—  , 
1  s  S-0 1  J  s-02  ' 


where 


(a1+w1+Si)X 

X+p+Xp 


s101+s1(l+X)+w1  (ai+Wi+sx)X 

2p  +  2p0^ 


and  H3  = 


-Si02  -  sx(l+X)  -  wx  (ai+w1+s1)X 
2p  2pQ~2 


The  transform  of  mw(t)  has  the  form: 


Fw  =  <s+^)Fs  -  si’ 


Direct  calculation  shows  that  H!+H2+H3  =  St  giving  the  following 


form  for  F.  : 
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w 


=  H.M  I  * 


+  H3(62+M^rQ: 


The  transform  of  m,(t)  has  the  form: 

a 

3  1 

F  =  — r—  +  F„,  which  can  be  put  into  the  form: 

a  S+A  S+A  s 

Him  1  h2m  ,  h3M  i 

r  -  _  £  +  -  - - +  -  — - — 

a  A  s  A+Qj  s-Qj  A+92  s-02 

Him  H2m  H3m  1 
+  ”  ~A~  X+07  ’  *+02  ^  s+A  ‘ 

It  can  be  shown  that  the  coefficient  of  is  equal  to  zero. 
Inversion  of  the  transforms  given  above  yields  the  Case  3  solutions 
in  the  following  form. 


ms(t)  =  Hi  +  H2e9lt  +  H3e02t 

mw(t)  =  mHi  +  H2(01+M)e0lt  +  H3(02+M)e02t 


ma(t) 


HiM  H2m  r  H3m  t 

-  +  -  e0il  +  -  e02t 

A  A+0i  A+02 


The  results  given  in  this  Chapter  are  rather  complicated,  and  it  is 
desirable  to  check  their  correctness  before  using  them  in  numerical 
analyses  such  as  those  given  in  Chapter  4.  Two  types  of  tests  on 
solution  correctness  have  been  made.  For  each  case  it  has  been  shown 
by  direct  calculation  that  m  (t)+m  (t)+m  (t)  =  a^Wx+Sj  as  required  by 

d  W  ♦  b 

the  conservation  assumption.  A  more  sophisticated  test  of  solution 
correctness  is  obtained  by  noting  that  both  Case  1  and  Case  3  solutions 
should  tend  to  the  Case  2  solutions  as  p-»0.  By  allowing  (5  to  approach 
zero  in  the  Case  1  and  Case  3  solutions,  it  is  found  that  the  Case  2 
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solutions  are  indeed  obtained.  Though  the  limit  calculations  are  detailed 
to  carry  out,  they  are  direct  calculations  and  are  omitted  from  this 
report. 

It  should  be  noted  that  the  initial  value  of  time  t  has  been  taken 
as  t  =0  in  deriving  the  above  forms  of  solution.  This  means  that  in 
model  analyses,  the  time  variable  must  be  redefined  each  time  the  model 
changes  region  form.  This  is  discussed  more  fully  in  Chapter  4. 

CHAPTER  3.  REGION  2  SOLUTIONS 


As  specified  in  Section  1.2,  the  Region  2  form  of  the  bulk  availa¬ 
bility  model  is  given  by  the  value  K(t)  =  r  +  (p-l)m (t).  This  results 

b 

in  the  following  system  of  differential  equations: 
ma(t)  =  pms(t)  -  Ama(t), 
mw(t)  =  Ama(t)  -  r  -  (p-l)ms(t), 
mg(t)  =  r  -  mg(t) . 

Initial  conditions  are  specified  for  some  value  of  time,  t  ,  at  which 
the  model  behavior  becomes  governed  by  the  Region  2  form.  As  in 
Chapter  2,  the  initial  values  are  defined  to  be  ma(tQ)  =  alt  ms(to)  =  si» 
and  mw(t0)  =  wx.  Because  of  the  conservation  assumption,  the  Region  2 
solutions  also  satisfy  the  condition: 

ma(t)  +  mg(t)  +  mw(t)  =  ax  +  sx  +  wx  for  all  t. 

The  equation  for  m  (t)  is  not  coupled  to  the  quantities  m  (t)  and 

b  Q 

mw(t)  so  that  it  can  be  solved  directly  for  mg(t).  It  is  a  first  order 
linear  differential  equation  whose  solution  is: 
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m  (t)  =  +  r(l-e_t). 

s 


This  value  can  be  put  into  the  equation  for  m_(t)  to  obtain  the 

a 


equation : 


ma(t)  =  pr  +  (j(s!-r)e  -  Ama(t). 


This  equation  is  solved  directly  to  obtain: 


,<t)  =  £  a-e-«)  ♦  (e-1  -  e-M)  ♦ 


Now  using  the  values  already  obtained  for  ms(t)  and  ma(t),  the 
equation  for  mw(t)  becomes: 

(jr-Apst  _xt  (Si-r)  _ 

mw(t)  =  [-j^i -  +  ^ai]e  +  A_1—  (A+m-1)  e  *•. 


Direct  integration  of  this  differential  equation  yields: 

ur  i  pr-Apsj  _x 

mw(t)  =  a^Wi+Si  -  -  r  -  j-  [  A_r-  +  Aaje 


(si-r)(A+p-l)  _t 
A-l  e  • 


As  a  check  on  these  results,  direct  calculation  shows  that  the 

conservation  condition  holds  for  the  Region  2  solutions  given  above. 

It  should  be  noted  that  the  initial  value  of  time  t  has  been  taken 

o 

as  t  =0  in  deriving  the  above  forms  of  solution.  This  means  that  in 
model  analyses,  the  time  variable  must  be  redefined  each  time  the  model 
changes  region  form.  This  is  discussed  more  fully  in  Chapter  4. 


Because  of  the  simple  form  in  which  the  Region  2  solutions  depend 
on  t,  the  limit  values  as  t  increases  are  easily  obtained.  They  are: 
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limit  m  (t)  =  r, 

t-*»  s 

limit  m,(t)  =  4^ 

t-x» 

limit  m,  ,(t)  =  ai+Wi+Si  -  ^  -  r. 

t-*» 

It  should  be  observed  however  that  as  t  increases,  the  model  may 
not  remain  in  Region  2.  If  it  does  not,  then  the  limit  values  are  unde¬ 
fined.  The  condition  for  the  model  to  be  in  Region  2  is  that  mw(t)  £ 
r  +  ((j-l)m  (t).  The  limit  form  of  this  condition  requires  that  ax+Wi+Si 
z  pr  +  ^  +  r. 

This  is  a  condition  that  relates  all  six  parameters  of  the  model  to 
give  a  set  of  valups  that  result  in  the  Region  2  limit  forms  for  the 
model  descriptors  m  ,  mg,  and  n^. 

Though  numerical  evaluation  is  required  for  typical  model  analyses 
in  either  Region  1  or  Region  2,  some  special  cases  can  be  considered 
directly.  As  an  illustration  of  such  a  consideration,  the  following 
special  case  has  been  studied. 

Consider  the  Region  2  solution  in  which  s,  =  r  so  that  the  service 
facility  starts  fully  utilized.  Then  m  (t)  =  r  and  the  facility  continues 
to  be  fully  utilized.  In  this  case,  the  remaining  units  are  distributed 
between  the  available  and  waiting  states  according  to  the  functions: 

V*)  *  “x 

^w(t)  =  ax  +  Wj  -  ^  -  (ax  -  )  e'At. 


The  condition  for  this  to  govern  the  distribution  of  units  is  that 
mw(t)  t  Mr;  if  this  condition  is  not  satisfied,  the  system  will  move  into 
the  Region  1  solution.  In  particular,  for  a  steady  state  limit  solution  to 
exist  in  this  form,  it  is  required  that  a^Wj  ^  pr(l+^).  It  is  an  indi¬ 
cation  of  the  intrinsic  non-linearity  of  the  model  that  the  steady  state 
depends  on  initial  values.  An  example  illustrating  the  situation  when 
the  condition  is  satisfied  is  given  in  Section  4.2. 

CHAPTER  4.  NUMERICAL  EXAMPLES  OF  SOLUTIONS 


The  bulk  availability  model  has  been  formulated  so  that  its  des¬ 
criptive  functions  m,(t),  m0(t),  and  m,  ,(t)  have  different  forms  in  the 
solution  Regions  1  and  2.  As  the  system  moves  from  one  region  of 
solution  to  another,  time  must  be  reinitialized  and  the  initial  distribution 
of  units  among  the  three  states  of  the  model  must  be  specified  by  the 
values  obtained  in  the  previous  region.  The  solutions  themselves  have 
been  obtained  in  explicit  form  but  are  rather  complicated.  Certainly 
they  are  too  involved  for  hand  calculation.  In  this  chapter,  a  computer 
implementation  is  described,  and  it  is  used  to  obtain  some  example 
results.  These  examples  illustrate  the  use  of  the  model  for  bulk  avail¬ 
ability  analyses  and  show  the  kind  of  results  obtained  from  the  mathe¬ 
matical  solutions.  A  major  feature  of  interest  in  model  analysis  is  the 
relation  of  the  system  descriptors  rn  (t),  nri  (t),  and  m  (t)  to  the 
system  parameters  p,  r,  and  model  parameters  alf  slf  and  wx . 

Section  4.1  Computerized  Solution  Calculation 

As  pointed  out  in  each  Region  solution  section,  the  solution  forms 
given  are  based  on  an  assumption  that  the  initial  time  value  is  zero.  In 


Page  26 


a  real  solution  where  the  model  form  may  pass  from  one  solution  form  to 
another,  it  is  necessary  to  use  two  time  values.  One  value  of  time 
represents  the  actual  time  while  the  other  time  value  is  the  time  within 
a  particular  region  solution.  Each  change  from  one  region  form  to 
another  requires  an  initialization  of  the  "local"  time  variable  to  zero. 
When  this  is  done,  the  actual  time  variable  is  incremented  by  a  value 
equal  to  the  local  time  value  just  prior  to  its  initialization  and  the 
change  to  a  new  region  form  solution.  Depending  on  the  six  parameters 
that  define  the  model,  three  types  of  solutions  can  arise: 

(i)  the  solution  can  start  in  one  region  and  remain  in  that  region 
for  the  duration  of  the  analysis; 

(ii)  the  solution  can  start  in  one  region  and  move  to  the  other 
which  then  determines  its  value  for  the  duration  of  the  analysis; 

(iii)  the  solution  can  move  from  one  region  to  another  then  back 
again,  continuing  to  change  regions. 

A  combination  of  (ii)  and  (iii)  can  also  occur,  but  the  major  types 
of  behavior  are  given  by  these  forms  of  solution.  Thus  a  computerized 
model  analysis  must  be  formulated  in  such  a  way  that  it  can  represent 
each  of  these  distinct  kinds  of  system  behavior. 

In  addition  to  keeping  track  of  actual  time,  it  is  necessary  to 
match  the  values  of  m  ft),  rn  (t),  and  m,  ,(t)  at  the  end  of  one  model 
form  situation  with  the  initial  values  of  these  quantities  at  the  next 
model  form  when  a  change  of  form  is  required. 

Figure  5  shows  a  flow  chart  of  how  the  computerized  model 
analysis  is  carried  out.  To  keep  track  of  time,  three  different  "time 
like"  variables  are  used.  The  variable  T  measures  time  from  the  start 
of  one  Region  form.  It  changes  in  increments  of  specified  amount  DT. 
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The  value  of  DT  must  be  selected  rather  carefully  because  it  governs 
the  time  values  at  which  the  solution  forms  are  evaluated.  If  the 
solution  remains  in  a  single  Region,  DT  can  be  relatively  large,  its 
value  subject  only  to  the  detail  of  solution  form  required  in  a  particular 
analysis.  However,  when  the  solution  changes  Regions,  it  should  do  so 
as  near  to  the  critical  values  of  the  solution  as  possible.  If  DT  is  too 
large,  the  solution  from  one  Region  will  continue  into  the  other  Region 
on  transition  and  the  "initial  values"  will  be  incorrect  being  based  on 
the  previous  solution  form  which  has  ceased  to  be  appropriate  at  some 
previous  time.  This  effect  can  introduce  significant  errors  in  a  model 
analysis.  The  most  appropriate  values  for  DT  in  a  specific  analysis  can 
be  determined  by  a  try-and-test  procedure.  Fortunately,  a  computer 
run  of  the  model  does  not  require  extensive  time  or  storage  so  it  is 
feasible  to  try  several  values  of  DT  to  insure  a  realistic  set  of  results. 
The  variable  TOTT  records  each  of  the  times  spent  by  the  system  in 
each  solution  Region.  Each  time  there  is  a  change  of  region,  the  value 
of  TOTT  up  to  the  start  of  that  Region  solution  is  increased  by  the 
final  value  of  T  for  that  solution.  Actual  time  is  given  by  the  variable 
TLOC  which  is  the  sum  of  the  present  value  of  TOTT  and  the  value  of 
T  obtained  so  far  in  the  present  Region  solution  at  the  termination  of 
the  analysis.  An  analysis  is  terminated  when  TLOC  reaches  a  value  TF 
specified  by  the  program.  The  quantities  DT  and  TF  are  input  values 
that  govern  the  detailed  calculation  of  the  solutions,  they  are  program¬ 
ming  parameters  rather  than  system  parameters. 

The  solution  Region  that  should  be  used  is  determined  by  the 
condition  m  ( t )  <  r-m„(t)+|jrn  (t)  as  discussed  in  Section  1.2 

w  ,  s  S 

(Figure  2).  The  proper  Region  is  designated  by  a  flag  variable  KR 
which  takes  the  values  1  and  2  to  denote  Regions  1  and  2  respectively. 
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At  each  iteration  of  the  solution  evaluation  procedure,  the  con¬ 
dition  on  mw(t)  is  checked  to  see  if  a  change  of  Region  should  take 
place.  When  a  change  of  Region  takes  place  TOTT  is  updated,  the 
values  of  in  (T),  rn  (T),  and  m,  ,(T)  are  used  as  the  intial  values  for 
the  next  Region  solution  form,  and  T  is  initialized  to  DT  as  a  starting 
value  for  the  new  Region  solution.  This  procedure  is  indicated  in 
Figure  5.  In  the  detailed  program,  a  listing  of  which  is  given  in  the 
Appendix,  care  must  be  taken  in  evaluating  the  Region  1  solutions.  In 
that  region,  the  three  Cases  give  different  forms  of  solutions  as  des¬ 
cribed  in  Chapter  2.  However,  the  notation  defining  the  quantity  0 
must  be  properly  interpreted  in  each  case.  This  is  particularly  true  in 
Case  3  where  D  <  0  so  that  02  is  in  fact  negative  making  a2-02  a 
positive  quantity  for  all  values  of  a  and  0.  The  detailed  program  must 
test  for  each  case  and  interpret  the  numerical  signs  so  as  to  give  the 
correct  results.  It  must  also  determine  how  often  solution  values  are  to 
be  given  as  output  and  what  detail  of  output  information  is  required. 
These  features  are  indicated  by  additional  program  parameters  as 
described  in  the  Appendix. 

In  the  next  section,  computerized  analyses  of  several  examples  are 
given.  They  were  obtained  using  the  program  given  in  the  Appendix 
which  follows  the  flow  diagram  of  Figure  5.  Examples  such  as  these 
require  very  short  running  times  and  little  storage  for  their  execution. 

Section  4.2  Examples 

To  illustrate  the  kinds  of  solutions  that  can  be  obtained  from  the 
model,  some  examples  are  given  in  this  section.  First,  a  simple  example 
is  given  for  which  the  solution  form  can  be  expressed  explicitly.  Then 


Input: 

DT.TF 
Al,  Wl,  SI 
M.  A,  R 


Figure  5 

Flow  Diagram  of  Computational  Program 
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some  numerical  examples  are  given  that  employ  the  computer  evaluation 
procedure  described  above  in  Section  4.1. 

Example  1.  Let  A  =  1,  p  =  2,  r  =  10,  and  assume  initial  values 
ax  =  70,  sx  =  10,  and  wx  =  20.  The  constant  value  of  total  population 
of  units  is  therefore  100.  The  condition:  wx  <  r-s1+ps1  becomes 
20  <  10-10+20  which  is  not  satisfied  so  a  Region  2  type  of  solution  is 
obtained.  In  Chapter  3,  the  limit  behavior  of  the  model  when  sx  =  r 
was  discussed  and  it  was  found  that  the  solution  would  remain  in 
Region  2  provided  the  condition  a!+wx  £  pr(l+^)  was  satisfied.  In  the 
present  example,  this  condition  has  the  values  70+20  £  10(1+1)  and  is 
satisfied.  Therefore,  the  solution  starts  in  Region  2  and  remains  in 
Region  2.  The  explicit  solution  is: 

ma(t)  =  20  +  50e-t 

mw(t)  =  70  -  50e_t 

ms(t)  =  10 

The  total  number  of  units  is  maintained  at  100  and,  in  the  steady  state, 
20  are  available,  10  are  in  service,  and  70  are  waiting  for  service. 

A  modification  of  the  initial  conditions  for  this  example  yields  quite 
different  results.  Let  ax  =  80,  sx  =  10,  and  wx  =  10.  Then  the  con¬ 
dition  wx  <  r-Sj+pSi  is  satisfied  in  the  form  10  <  20  and  the  solution  is 
in  Region  1  with  D  =  -4,  p2  =  -1,  and  a  =  2.  Thus  Case  1  type  solu¬ 
tions  must  be  used  to  represent  the  descriptive  functions,  which  have 
the  explicit  form: 

m$(t)  =  20  -  lOe"^1"  cost  -  30e”^t  sint 

ma(t)  =  40  +  40e~^  cost  +  20e"^1'  sint. 

mw(t)  =  40  -  30e"^  cost  +  10e~^  sint. 


The  conservation  condition  ma(t)+ms(t)+mw(t)  =  100  holds,  as  it 
should,  for  tnese  solutions. 

The  Case  1  examples  are  interesting  because  the  D  <  0  condition 
permits  oscillations  to  occur  causing  the  level  of  availability,  for 
example,  to  move  up  and  down  about  a  steady  state  value.  This  is 
illustrated  by  putting  the  availability  function  m=(t)  into  the  following 

d 

form: 

nr  (t)  =  40  +  20  JS  e”^  sin(t+0) 

a 

2 

where  sin  0  =  — .  In  numerical  terms,  the  function  is  expressed  as 
V5 

m.(t)  =  40  +  44.72e~^  sin(t+l.l). 

a 

The  limit  value  is  40  and  the  values  of  ma(t)  undergo  damped 
oscillation  about  this  value  as  shown  in  Figure  6.  The  limit  values  of 
mg(t)  and  mw(t)  are  seen  to  be  equal  to  20  and  40  respectively. 

The  remaining  examples  of  this  section  are  given  in  graphical  form 
where  the  function  values  are  obtained  by  using  the  computer  proce¬ 
dure  indicated  in  Figure  5.  Numerical  values  for  each  example  are 
given  in  the  Appendix. 

Example  2.  To  illustrate  a  Region  1  solution,  let  (j  =  .4,  \  =  .01, 
r  =  10,  and  set  the  initial  values  equal  to  al  =  89,  sx  =  5,  w1  =  6so 
that  the  total  number  of  units  is  100.  For  this  example,  D  =  .33  >  0  so 
it  has  a  Case  3  solution.  The  values  of  ma(t),  ms(t),  and  mw(t)  are 
shown  in  Figure  7  for  20  time  units.  Because  the  solution  does  not 
change  regions,  the  time  increment  used  in  the  evaluation  program  can 
be  large.  A  value  of  DT  =  1  was  used  for  the  numerical  values  shown 
in  Figure  7. 
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Bulk  availability  of  the  kind  shown  in  this  example  can  be  charac¬ 
terized  as  a  rather  pure  or  classical  type  of  availability  behavior  as 
distinct  from  a  classical  reliability  or  mixed  forms  illustrated  in  other 
examples.  The  service  rate  is  greater  than  the  failure  rate  by  a  large 
enough  amount  that  the  availability  function  is  able  to  reach  a  steady 
state  value  of  approximately  97  units.  Of  the  remaining  three  units, 
two  are  in  service  and  one  is  waiting  for  service,  in  the  steady  state. 

Example  3.  A  Region  2  solution  is  illustrated  in  this  example. 
The  parameter  values  are:  p  =  .2,  k  =  .2,  and  r  =  10.  Initial  values 
are  taken  to  be  a:  =  70,  sx  =  10,  and  wT  =  20  giving  a  total  of  100 
units.  Here,  the  initial  number  in  service,  slt  is  equal  to  the  capacity 
of  the  service  facility.  Under  these  conditions  with  equal  values  for 
service  and  failure  rates,  the  larger  number  of  units  initially  available 
forces  the  size  of  the  population  of  units  waiting  for  service  to  increase. 
This  continues  until  a  steady  state  is  reached  with  in  (t)  =  in  (t)  =  10 
units  and  with  80  units  in  the  waiting  state.  The  values  are  shown  in 
Figure  8.  Since  the  solution  remains  in  Region  2,  a  relatively  large 
value  of  DT  can  be  used.  In  this  example  DT  =  1.  The  model  acts 
like  a  classical  reliability  type  model  in  that  the  failure  rate  is  strong 
enough,  relative  to  the  capability  for  service,  that  the  lack  of  available 
units  falls  until  a  relatively  low  steady  state  level  is  reached. 

Example  4.  This  example  illustrates  an  extremely  rapid  transition 
from  a  Region  1  to  a  Region  2  type  of  solution.  The  parameter  values 
are:  p  =  .4,  \  =  .4,  and  r  =  10.  The  initial  values  are:  ax  =  89, 
sa  =  5,  and  wx  =  6  for  a  total  of  100  units.  Because  of  the  rapid 
transition  between  regions  a  small  value  of  DT  is  required  so  as  to  not 
miss  or  greatly  distort  the  transition  affect.  A  value  of  DT  =  .001  time 
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units  was  used  for  this  example.  The  values  shown  in  Figure  9  give 
the  Region  1,  Case  1  (D  =  -.6<0)  solution  which  shifts  to  a  Region  2 
solution  at  T  =  .035  time  units.  The  values  of  the  descriptive  functions 
m„  (t),  m  ft),  and  mIT (t)  do  not  change  much  over  this  short  time 
interval.  Though  the  example  is  given  to  illustrate  how  fast  the  model 
can  shift  from  one  solution  form  to  another,  it  should  be  noted  that  in 
applications,  the  time  units  may  be  such  as  to  give  practical  interest  to 
relatively  short  time  periods.  This  is  discussed  more  fully  in  Chapter  5. 
The  most  important  aspect  of  this  example  is  the  need  for  using  a  very 
small  DT  value.  Otherwise  the  true  nature  of  the  solution  is  lost  and  a 
numerical  solution  is  obtained  that  does  not  properly  represent  the 

nature  of  the  model  solution  for  this  set  of  parameters  and  initial 

conditions . 

Example  5.  A  less  rapid  transition  from  a  Region  1  to  a  Region  2 
solution  is  illustrated  by  this  example.  The  parameter  values  are 
p  =  .4,  \  =  .2,  and  r  =  10.  Initial  conditions  are  ax  =  100,  sx  =  0,  and 
Wi  =  0  so  that  at  the  start  of  the  analysis  all  units  are  available. 
Though  the  transition  between  regions  is  not  as  rapid  as  for  Example  4 
it  is  sufficiently  fast  that  care  is  required  in  the  selection  of  the  evalu¬ 
ation  time  interval  if  correct  Region  2  solutions  are  to  be  obtained. 
The  value  used  is  DT  =  .05  time  units.  Transition  from  Region  1  to 
Region  2  takes  place  at  T  =  .5  time  units  and  by  2  time  units,  the 

service  facility  has  almost  reached  its  capacity  of  10  units.  Though  the 
example  does  not  go  completely  to  a  steady  state  in  2  time  units,  the 
model  descriptors  have  come  close  to  steady  state  values  by  this  time. 
Figure  10  shows  the  results  and  indicates  extrapolated  steady  state 

values  of  ma(t)  =  68,  mg  =  10,  and  mw  =  22. 


Region 


Time  Units 
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CHAPTER  5.  BULK  AVAILABILITY  WARRANTIES 

Over  the  past  several  years,  customers,  such  as  the  U.S. 
Government,  have  become  increasingly  concerned  with  the  satisfactory 
operation  of  equipment  after  it  has  been  delivered.  Acceptance  test 
procedures  insure  that  equipment  meets  contractual  requirements  at  the 
time  of  delivery  but  provide  little  indication  of  successful  subsequent 
operation.  Reliability  tests  that  form  part  of  many  such  procedures 
have  minimal  predictive  value  due  to  such  factors  as  lack  of  statistically 
significant  data,  intrinsic  variability  in  operational  conditions,  and 
limited  periods  for  observation  of  equipment  before  completion  of  accep¬ 
tance  procedures.  This  situation  has  led  to  increased  interest  in 
warranty  policies  that  may  be  included  as  part  of  contract  specification. 
In  particular,  the  concept  of  the  Reliability  Improvement  Warranty 
(RIW)  provides  a  means  for  the  contractor  to  assume  some  of  the  risk 
associated  with  the  customer's  acceptance  of  equipment  which  subse¬ 
quently  fails  to  operate  properly.  The  development  of  RIW  models, 
their  analysis,  and  their  application  to  contract  definition  has  brought 
about  a  well-defined  field  of  study  blending  classical  warranty  concepts 
with  stochastic  and  other  mathematical  techniques  from  the  field  of 
Reliability  (7). 

Such  warranties  may  be  considered  from  several  points  of  view. 
They  may  be  considered  as  motivation  for  a  contractor  to  strive  for 
high  reliability  products.  They  can  be  a  means  for  clearly  specifying 
the  division  of  risk  assumption  between  contractor  and  customer  due  to 
product  failure.  They  can  aslo  provide  a  range  of  design  cost  trade¬ 
offs  in  which  a  contractor  can  choose  between  spending  money  on 
producing  a  more  reliable  product  or  on  paying  warranty  costs  (4). 
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An  immediate  extension  of  the  RIW  concept  is  to  consider  war¬ 
ranties  from  the  point  of  view  of  availability  rather  than  reliability. 
This  places  the  focus  on  the  availability  of  a  product  for  use  when 
needed  which  can  be  achieved  by  improving  its  reliability  and/or 
providing  improved  repair/replace  options.  Availability  warranties  seem 
particularly  attractive  to  the  customer.  However,  they  may  be  viewed 
as  undesirable  contract  obligations  by  a  contractor.  By  combining  the 
concepts  of  incentive  contracts  and  availability  warranties  a  wide  range 
of  contract  structures  can  be  formulated  in  which  a  combination  of  risk 
assumption  and  profitability  can  be  achieved  (8). 

Much  of  the  study  of  availability  warranties  has  considered  the 
availability  of  a  single  item  or  system  which  is  subject  to  failure  and 
repair/replace  service  over  the  duration  of  the  warranty  period  (9), 
(3).  This  represents  a  classical  viewpoint  in  the  study  of  availability. 
The  present  paper  deals  with  another  kind  of  availability  consideration 
which  considers  a  large  number  of  units  rather  than  a  single  item. 
Each  of  the  single  units  is  relatively  minor  in  its  required  function  and 
availability  is  specified  by  having  a  minimum  number  of  units  available. 

In  this  study,  the  availability  concepts  are  those  of  bulk  avail¬ 
ability  rather  than  individual  item  availability.  Therefore  a  determi¬ 
nistic  model  for  bulk  availability,  as  defined  in  Chapter  1,  has  been 
used  as  a  basis  for  anlaysis  rather  than  the  more  common  approach 
based  on  stochastic  models.  This  procedure  seems  analogous  to  the  role 
of  thermodynamics  in  describing  gross  characteristics  of  matter  that 
require  statistical  mechanics  for  their  detailed  investigation.  This 
analogy  provided  the  guidelines  for  development  of  the  present  bulk 
availability  model. 
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Bulk  availability  deals  with  such  things  as  a  large  work  force 
(manpower  availability),  fleets  of  equipment,  or  graceful  degradation  of 
complex  systems.  An  example  is  a  military  electronic  communications 
device  which  may  be  supplied  in  several  thousand  units.  A  particular 
group  may  have  100  such  devices  assigned  to  it  and  require  80  for 
completely  satisfactory  operation.  This  would  allow  20  of  the  devices  to 
be  either  undergoing  service  or  waiting  for  service.  If  more  than  20 
were  out  of  order,  the  bulk  system  would  be,  to  some  degree,  unavail¬ 
able.  If  an  availability  warranty  is  in  effect  in  such  a  situation,  the 
contractor  would  have  to  pay  penalty  costs  specified  by  the  warranty 
contract.  Such  costs  can  be  viewed  in  various  ways:  as  out-and-out 
payment  of  penalty  dollars,  as  a  requirement  to  spend  money  on 
improved  repair/replace  facilities,  or  as  the  free  (or  reduced  cost) 
supply  of  one  or  more  additional  units  by  the  contractor. 

This  report  has  described  the  deterministic  bulk  availability  model 
and  how  analyses  can  be  carried  out  using  the  model.  It  gave  some 
examples  of  the  behavior  of  the  bulk  availability  models.  How  they 
specify  availability  as  a  particular  level  of  available  units  will  be 
described  in  this  chapter. 

Bulk  availability  warranties  are  defined  as  costs  incurred  by  the 
contractor  when  the  number  of  available  units  falls  below  contractually 
specified  levels.  The  cost  trade-offs  available  to  a  contractor  are 
illustrated  by  the  examples.  Conclusions  about  the  utilization  and 
applicability  of  the  bulk  availability  warranty  concept  are  discussed  in 
the  final  chapter  of  the  report. 

The  most  immediate  application  of  the  model  to  a  bulk  availability 
warranty  analysis  is  when  a  penalty  fee  is  charged  to  the  contractor 
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whenever  the  number  of  available  units  falls  below  a  specified  level. 
This  kind  of  response  to  a  lowered  level  of  availability  has  no  effect  on 
the  time  functions  of  the  model.  When  necessary,  the  fee  is  paid  and 
the  model  continues  to  represent  the  bulk  availability  situation.  Pay¬ 
ments  of  this  kind  can  be  as  a  single  penalty  or  can  be  tied  to  the 
length  of  time  over  which  the  availability  is  too  low.  In  situations  such 
as  those  shown  in  Figures  8  and  10,  the  availability  function  decreases 
monotonically .  In  such  a  case  if  the  function  falls  below  the  specified 
availability  level  before  the  warranty  time  has  expired,  a  penalty  will  be 
incurred  by  the  contractor.  If  it  is  of  the  time  duration  type  it  is 
clearly  desirable  for  the  low  level  of  availability  to  be  reached  as  near 
to  the  end  of  the  warranty  period  as  possible.  This  shows  the  value  of 
a  time  dependent  penalty  to  the  customer.  It  can  also  be  used  tc  help 
the  contractor  in  both  product  design  and  contract  structure  decisions. 
In  a  situation  like  that  shown  in  Figure  7  the  availability  is  increasing. 
In  such  a  case,  the  contractor  might  incur  a  warranty  penalty  initially 
but  would  subsequently  be  free  of  warranty  payments.  This  indicates 
the  importance  of  design  considerations.  A  high  enough  initial  state  for 
availability  would  prevent  any  warranty  costs.  However,  the  cost  to 
achieve  such  an  initial  level  could  exceed  the  initial  warranty  cost  and 
therefore  be  undesirable.  Another  consideration  is  the  increasing 
availability  function.  This  kind  of  product  performance  may  be  better 
than  required,  at  a  cost  to  the  contractor.  An  alternative  design,  with 
lower  cost  and  lower  availability  might  serve  both  parties  better,  giving 
satisfactory  performance  to  the  customer  and  higher  profit  to  the  con¬ 
tractor.  Figure  6  shows  a  variation  in  the  availability  function.  In  a 
case  such  as  this,  the  proper  level  for  availability  may  be  in  effect  at 


some  periods  of  time  and  not  during  other  periods.  This  is  an  inter¬ 
esting  case  from  the  warranty  point  of  view.  It  allows  a  wide  range  of 
contract  stipulations  that  might  not  be  considered  without  a  knowledge 
of  expected  behavior  of  this  kind.  Thus,  in  this  type  of  bulk  availa¬ 
bility  model,  it  is  particularly  important  to  have  some  kind  of  qualitative 
guidelines  for  the  behavior  of  the  availability  function,  as  provided  by 
the  model  analysis. 

A  more  realistic  type  of  warranty  obligation  is  to  provide  some 
form  of  correction  for  the  occurrence  of  low  value  of  the  availability 
function.  When  the  availability  function  falls  below  a  specified  value 
the  contractor  may  be  required  to  bring  the  number  of  available  units 
up  to  a  level  that  will  yield  a  satisfactory  value.  This  will  require 
replacement  and/or  repair  of  units.  An  action  of  this  kind  interrupts 
the  time  functions  of  the  model  analysis.  When  it  is  desired  to  employ 
this  kind  of  response  to  a  warranty  obligation,  the  mathematical  model 
must  be  interrupted.  A  new  set  of  initial  conditions  are  defined  as  a 
combination  of  what  was  being  specified  by  the  model  upon  interruption 
and  actions  taken  by  the  contractor  before  resumption  of  the  model 
analysis.  Some  way  to  measure  time  must  also  be  selected  so  as  to 
represent  time  in  a  useful  way.  Most  likely  the  time  should  be  con¬ 
tinued  to  be  measured  from  the  beginning  of  the  warranty  period.  In 
this  case,  the  time  between  model  interruption  and  the  new  start  should 
be  included  in  the  total  time  description.  This  time  measures  the  time 
required  to  bring  the  availability  up  to  the  desired  level.  It  is  seen 
that  the  repair/replace  procedure  for  satisfying  warranty  obligations 
can  be  studied  within  the  framework  of  the  bulk  availability  model  of 
this  paper  by  making  the  kinds  of  modifications  to  the  analysis 
described  above. 
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CHAPTER  6.  CONCLUSIONS  AND  EXTENSIONS 


A  deterministic  model  for  bulk  availability  has  been  described  in 
this  report.  It  has  a  fixed  total  population  and  what  may  be  called  an 
r-server  type  repair  queue  system.  In  some  ways  the  model  appears 
similar  to  a  stochastic  queuing  model  with  these  characteristics.  How¬ 
ever,  the  bulk  model  has  gross  state  transition  values  rather  than 
individual  transition  effects.  This  makes  it  more  like  a  closed  system 
dam  type  process  than  a  particle  queue  process.  Though  the  availa¬ 
bility  level  curve  as  a  function  of  time  looks  much  like  a  mean 
availability  curve  for  various  queue  processes,  it  does  not  seem  to 
follow  in  any  direct  way  from  such  processes.  This  is  most  likely  due 
to  the  bulk  transition  effect  which  prevents  the  application  of  the  basic 
non-cascade  assumption  necessary  for  the  formulation  of  most  particle 
queue  models.  Research  in  a  dam  type  process  expressing  properties 
similar  to  the  present  model  might  be  of  theoretical  interest  but  is 
considered  outside  the  scope  of  the  research  described  here. 

The  model  description  and  analysis  described  in  this  report  shows 

that: 

•  The  bulk  availability  model  can  be  employed  as  a  useful  tool  in 
trade-off  analyses  between  manufacturing,  logistic,  and  warranty  costs. 
By  using  several  versions  of  a  proposed  model,  with  different  parameter 
values,  a  set  of  possible  availability  functions  can  be  generated.  This 
kind  of  study  can  give  guidelines  for  determining  the  best  values  or 
ranges  of  values  to  aim  for  in  designing  a  product  to  meet  selected  cost 
trade-off  values. 

•  The  deterministic  model  is  simpler  to  formulate,  interpret,  and 
explain  than  corresponding  stochastic  models.  This  is  particularly  true 


in  the  present  case  where  the  appropriate  stochastic  model  seems  to  be 
of  the  finite  dam  type  in  which  the  mathematics  is  more  involved  than 
for  queue  type  models.  Ease  in  application  and  in  intuitive  relevence 
between  parameter  values  and  results,  such  as  availability  levels,  is 
particularly  useful  in  contract  negotiations  where  one  party  may  be 
required  to  explain  or  justify  its  analyses  as  part  of  its  justification  for 
cost  trade-off  decisions. 

•  The  closed  form  of  direct  solution  for  this  kind  of  model  allows 
the  consideration  of  numerical  results  over  any  desirable  time  period, 
subject  only  to  considerations  of  transition  between  solution  regions. 
The  examples  given  in  Chapter  4  illustrate  the  ease  with  which  numer¬ 
ical  solutions  can  be  obtained  using  readily  available  computing 
resources. 

Though  the  calculations  required  are  sufficiently  involved  that  the 
model  could  not  be  used  without  computer  implementation,  the  computer 
time  and  storage  requirements  are  small.  For  example,  20  steps  of  the 
model  run  in  what  appear  to  be  no  time  to  an  interactive  terminal  user 
on  a  standard  computer  installation. 

•  The  bulk  availability  model  can  be  given  broader  interpretations. 
For  example,  if  it  is  to  be  used  as  a  model  for  manpower  availability  the 
quantity  A.  does  not  represent  manufacturing  effort,  it  relates  to  the 
profile  of  people  recruited  into  the  manpower  pool.  Thus  A  can  be 
related  to  recruitment  effort  costs  such  as  extent  of  physical  and  mental 
examination  of  candidates.  It  can  also  be  related  to  the  level  of 
physical  or  mental  condition  that  will  be  used  in  accepting  or  rejecting 
candidates. 

The  bulk  availability  model  given  in  this  report  could  be  made 
more  widely  applicable  if  it  was  extended  to  allow  time  variation  in  some 
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or  all  of  the  model  parameters  \,  p,  and  r.  In  the  present  model,  the 
most  usual  form  of  solution  is  to  move  to  Region  2  and  remain  there 
tending  toward,  and  usually  reaching  rather  soon,  a  steady  state  condi¬ 
tion  with  the  service  facility  full.  If  the  parameters  change  with  time  a 
more  dynamic  model  would  result  which  would  allow  transitions  back  to 
Region  1.  Such  a  dynamic  model  is  felt  to  be  a  reasonable  one  in 
modeling  situations  in  which  a  contractor  can  change  parameter  values 
or  operating  conditions  cause  them  to  change  during  the  warranty 
period. 
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APPENDIX 


Computer  Program 

A  listing  of  the  FORTRAN  computer  program  used  to  carry  out  the 
numerical  evaluation  of  the  model  is  given  below.  It  is  one  particular 
realization  of  the  flow  diagram  shown  in  Figure  5.  The  main  reason  for 
giving  this  listing  is  to  illustrate  the  simple  programming  required  to 
implement  the  model  in  terms  of  numerical  calculations. 


WR ITE ( 6 , *  >  'INPUT  DT  AND  TF ' 

READ( 5 , * )  DT.TF 

WRITE ( 9 , 300 )  DT.TF 

WRITE! 6 , *  )  'INPUT  A1 ,W1 ,  AND  SI* 

READ! 5 , *  )  A1.W1.S1 

WRITE (9,310)  A1.W1.S1 

WRITE! 6 , *  )  'INPUT  XMU.XLM.AND  R' 

READ! 5 , *  )  XMU , XLM , R 
WRITE! 9,320)  XMU. XLM, R 
WRITE ( 6 , *  )  'INPUT  K' 

READ! 5 , *  )  K 
WR I TE ( 9 , 360 )  K 
T=0. 

TQTT=0. 

TLOC-0. 

KTR-0 

AA=R-S 1 +XMU*S 1 
IF(Wl.GT.AA)  GO  TO  20 
KR=  1 

60  CONTINUE 

D=( 1 .-XLM)*( 1 .-XLM)-2.*XMU*( 1 . +XLM >+XMU*XMU 
ALP-! 1 . +XLM+XMU  )/2 . 

AB=ABS( D/4  .  > 

BTA«SQRT(AB ) 

C  WRITE!  9,350)  D  ,ALP  ,AB  ,BTA 

IF! D.EQ.0. )  GO  TO  70 
IF! O.GT.0. )  GO  TO  80 
CO!I  =  (  ALP*ALP  >  +  (  BTA*BTA  ) 

A=( S1+W1+A1 ) /COM 

Cl”!  XMU‘*S  1  *(  1 .+XLM)-XLM*(A1+W1  )  >*EXP ( -ALP*T ) 

C2»XMU*S1*( 1 .+XLM*XLM-XMU-XMU*XLM)+W1',!2.*XMU-XLM*XLM 
1  +!!LM+XLM*XMU  ) 

C3«-XLM*A1*( 1  ,+XLM  +  XMU  ) 

XMS»XLM*A+(Cl*COS( BTA*T )/COM >+( ( CZ+C3  )*E XP( -ALP*T )*S I N! BTA*T> 
1/(2. *BTA*COM )  > 

C4*( -(  S 1 +W1 >*XMU+A1 *XLM*( 1 . +XMU  > )*EXP( -ALP*T  > 

C5*S1 *XMU* ( XLM+XMU+2.*XLM*XMU-1 .  >-Wl*XMU*( 1 . +XLM+XMU > 

C6«A1*XLM*! 1 ,+XMU*XMU-XLM-XLM*XMU > 

XMA-XMU*A+(C4*COS( BTA*T )/COM )+! ( C5+C6 )*EXP! -ALP*T )*S IN! BTA*T ) 
1/(2. *BTA*COM  >  ) 

C7»(W1*< XLM  +  XMU  >-XLM*XMU*( Sl+Al > > *EXP  ( -ALP *T > 
C8*W1*(XLM*XLM+XMU*XMU-XLM-XMU )+Al *XLM*<  2 . *XLM+XMU+XLM*XMU 
1 -XMU*XMU  ) 

C9*-S1*XLM*XMU*(  1 .  +XMU  +  XLM  > 

XMW«XMU»XLM*A-MC7*COS( 8TA*T l/COM )  +  ( ( C8+C9  >*EXP < -ALP *T >*S I N ( BTA*T > 
1/(2. *BTA*COM )  ) 

GO  TO  100 
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70  8  = !  S 1 +W 1 +A 1 )  / ! AL  P  *AL  P  ) 

C-W1+S1*! 1 . +XLM-ALP >-XLM*ALP*B 
Cl-CT'-'EXP!  -ALP*T) 

C2  =  !S1-XLM«8  >*£XP< -ALP*T> 

XMS*XLH'*B+C  1 +C2 

C3-! -S1-W1+XLM*! 2.*ALP-XLM)*B > 

XMA-XHU  *B  +  ( ( XMU*C3*EXP< -ALP*T >  >/< (XLM-ALP )*( XLM-ALP ) >  > 
l-<  XMU*C1/!ALP-XMU ) ) 

C4=!W1-XLM*XMU*B  >*EXP< -ALP*T ) 

Xny=XMU:‘XLM«B+C4  +  (  XMU-ALP  >*C1 
GO  TO  100 
30  TH1-BTA-ALP 
TH2--BTA-ALP 

H1-1A1+W1+S1  )*XLM/!  XLM+XMU+XLM*XMU  5 
C1=S1“< l . +XLM  )+Wl 

C2*(A1+W1+S1 >*XLM/( 2 .  *BTA  > 

H2-! < S1*TH1+C1 )/!2.*BTA))+!C2/THl > 

H3-!  <  -SI "TH2-C 1  >/<  2.*BTA>  >-!C2/TH2  > 
XMS-H1+H2*EXP!TH1*T)+H3*EXP!TH2*T> 

X;;W»XHU  ,'H1+H2*{TH1+XMU  >*EXP!  TH1*T)  +  H3*!TH2  +  XMU>*EXP!TH2*T> 
XMA"! H 1 -XMU/XLM  )  +  ( ! H2*XMU*E XP ( TH 1 *T > )/( XLM+TH 1 > ) 
l  +  ( ( H3*XHU*EXP( TH2*T ) >/( XLM+TH2 > ) 

GO  TO  100 
20  KR-2 
50  CONTINUE 

XMS"! S 1 -R  )*EXP! -T  )  +  R 
C 1 =XMU*<  S 1 -R )/! XLM-1  .  > 

XMA*( XMU*R*( 1 .-EXP< -XLM*T>  > /XLM >+C 1 *! EXP ! -T >-EXP ! -XLM*T >  > 
1+A1*EXP< -XLM"T) 

C2-A1+W1+S1-R-! XMU*R/XLM  > 

C3  =  ( XMU*! R-XLM*S 1 )/<  XLM-1 . )  +  ( XLM*A1  )  >/XLM 
C4-!S1-R>*!XLM+XMU-1. )/( XLM-1.  ) 

XNU-C2-C3*EXP< -XLM*T)-C4*EXP!-T) 

GO  TO  100 
100  TLOC-TUTT+T 
KTR-KTR+1 

IF(KTR.LE.K)  GO  TO  110 
WRITE (9,330)  TLOC , XMA , XMS , XMW 
WRITE (9,340)  KR.T.TOTT 
KTR-0 

110  CONTINUE 

IF(TLOC.GE.TF)  GO  TO  999 

IF( XMW. LT. R-XMS+XMU*XMS )  GO  TO  200 

I F ! KR . EQ . 2 )  GO  TO  150 

Al-XMA 

S l  -  XMS 

Wl-XMW 

TOTT-T+TOTT 

T-DT 

KR-2 

GO  TO  50 
150  T-T+DT 
GO  TO  50 

200  IF(KR.EQ.l)  GO  TO  250 
Al-XMA 
Sl-XMS 
Wl-XMW 
TOTT-T+TOTT 
T-DT 
KR- 1 

GO  TO  60 
250  T-T+DT 
GO  TO  60 
999  STOP 

300  FORMAT! 5X , ' DT- ' , F 10. 5 , 5X , ' TF - ' , F 10. 5 > 

310  FORMAT! 5X, 'Al-' , F10.5.5X, 'Wl-' , F10.5.5X, *S1"' , F10.5) 

320  FORMAT! 5X, * XMU- ' . F 10. 5 , 5X . ' XLM- ' . F 10. 5 . 5X . ’ R- ' , F 10 . 5  ) 

330  FORMAT! 5X. 'TLOC-' .F10.5.5X, 'XMA-' , F 10. 5 , 5X , ' XMS- ' ,F10.5, 
15X, ’XMW-' ,F10.5 ) 

340  FORMAT! 5X, ' KR- ' , 1 3 , 5X , ' T- 1 , F 10. 5 , 5X , 'TOTT-' ,FI0.5> 

C  350  FORMAT! 5X,‘ TEST  ',4F10.S> 

360  FORMAT! 5X, *K-  ’.IS) 

FNH 


Page  49 


Numerical  Data  for  Examples 

The  numerical  data,  produced  by  the  computer  program  listed 
above,  are  given  below  for  the  examples  given  in  Chapter  4  of  the 
report.  They  are  listed  by  Example. 


Example  shown  in  Figure  7: 


37=  1 .ore 7. 7  TF  =  2 0.00000 

Al*  83  .00000  W 1  *  £ . 00000  SI-  5 . 00000 

XMU  =  0.4  0C00  XLM=  0.01000  R*  13. 2  2  B  SOS 

TEST  0.33210  0.70500  0.0S303  0.2S814 

TLOC*  0.00000  XX  A*  39.00001  XMS-  5.00000 

KR*  1  T=  0.00000  TOTT  =  0.00050 

TEST  0.23210  0.70500  0.03303  0.2S314 

tlc:*  1  0  XMA*  90.56210  XMS*  6. 65261 

.  1  T=  1.00030  TOTT*  0.00000 

TEST  0.33210  0.70500  0.03303  0.2S814 

TLOC*  2.00000  XMA-  92.24326  XMS*  6.15701 

KR*  1  T  =  2.00000  TOTT  =  0.00000 

i'EST  0.32210  0.70500  0.03303  0.23814 

UCC-  3.00000  X MA  =  93.53267  XMS*  5.22056 

KR-  1  T*  3.00000  T0TT=  0.00000 

TEST  0.33210  0.70500  *0.08303  0.28814 

TLOC*  4.00000  XMA*  94.57150  XMS*  4.40021 

<R»  1  T  *  4.00000  TOTT*  0.00000 

TEST  0.33210  0.70500  0.03303  0.28314 

TLOC*  5.00000  XMA*  95.24976  XMS*  3.77152 

KR-  1  T  =  5.00300  TOTT*  0.00000 

TEST  0.33210  0.70500  0.03303  0.23314 

tloc*  6.00000  x:-;a=  95.70903  xms-  3.32699 

KR*  1  T*  6.00000  TOTT*  0.00000 

~t  ST  0.33210  0.70500  0.08303  0.23814 

TLOC*  7.00000  XMA-  96.01537  XMS*  3.02233 

KR*  1  T*  7.00000  TOTT*  0.00000 

TES”  0.33210  0.70500  0.03303  0.28314 

TLOC-  8.00000  XMA*  96.22060  XMS*  2.81322 

KR*  1  T*  8.00000  TOTT*  0.00000 

TEST  0.33210  0.70500  0.08303  0.28814 

TLOC-  9.00000  XMA*  96.35332  XMS*  2.68183 

KR*  1  T=  9.00000  TOTT-  0.00000 

TEST  0.33210  0.70500  0.03303  0.28814 

TLOC-  10.00000  XMA-  95.44513  XMS*  2.59135 

KR-  1  T*  10.00000  TOTT*  0.00000 

TE  ST  0.33210  0.70500  0.08303  0.23814 

TLOC*  11.00000  XMA*  96.50417  XMS*  2.53152 

KR*  1  T*  11.00000  TOTT-  0.00000 

TEST  0.33210  0.70500  0.0S303  0.28814 

TLOC*  1 2 .00000  XMA*  96.54308  XMS*  2.49200 

KR-  1  T-  12.00000  TOTT*  0.00000 

TEST  0.33210  0.70500  0.03303  0.28814 

TLCC*  13.00000  XMA*  36.53373  XMS*  2.45593 

KR-  1  T-  13.00000  TOTT-  0.00000 

TEST  0.33210  0.70500  0.03303  0.23814 

'-OC-  14.00000  XMa.-  96.58565  XMS*  2.44873 

KR-  1  7=  14.00000  TOTT*  0.00000 

TEST  0.332:0  0.70E00  0.03303  0.23814 

TLOC-  15.00000  XMA*  95.59579  XliS*  2.43739 

KR-  1  T »  15.00000  TOTT-  0.00000 

TEST  0.32210  0.70500  0.03303  0.28814 


XMW* 

6.00000 

XKW- 

2.77530 

XMW* 

1 .5SS75 

XMW* 

1.17677 

XMW* 

1.02830 

XMW- 

0.97873 

XMW- 

0.96393 

XMW* 

0.96030 

XMW* 

0.96119 

XMW* 

0.96236 

XMW* 

0.96346 

XMW* 

0.96432 

XMW* 

0.SS493 

XMW* 

0.96535 

XMW* 

0.96563 

XMW* 

0.96532 

.  ! 
V  1 


1 


2.42591 


XMV/»  2 . 56554 


TlCC* 
KR*  1 
TE  ST 

■  la. 2:222  XMA*  56.  5.741b  -  2.42591 

T  =  1 6  . i 62222  TOT  i  *  2 . !o 7.  ~  J 

2.232:2  2 . 7 2522  2.25322  .5.23314 

XMW* 

0.S65S4 

TLCC* 
KR=  l 
TEST 

17.22222  XMA*  95.55359  X:!S=  2.42499 

T  =  17.22222  TOTT*  2.22232 

2.22212  2.72520  2.0S203  2. 28314 

XMW* 

0.96602 

TLL'C* 
KR  =  1 

TEST 

19 .32202  XMA*  96.61219  XMS*  2.42174 

T  =  12.22222  TOT~-=  2.22202 

2.22212  2.72522  0.25322  5.23814 

XMW* 

0.96608 

TLOC* 
KR-  1 

13.00000  XMA*  S3. 51435  XMS*  2.41960 

T  =  19.00000  TOTT*  0 . 322:52 

XMW* 

0.96611 

Example  shown 

in  Figure  8: 

27*  I 

V  -?.“r  T/* 

•  4--viV..  ^ 

Tr*  20 

•  >-•  /.•  A-  J  ~ 

A  i  =  70 

'.22.220 

V/l-  20 

.22253 

SI*  10.' 

Y  /V  -r  f/ 

t  iS  At  A. 

XMU* 

0.22230 

XLM  = 

0.20323 

R*  10. 

.00000 

TLCC  = 

2. 22 2 22 

XMA* 

70.02220 

XM“  = 

10.00000 

XMW* 

20.00000 

KR*  2 

0.2023 Z 

TOTT* 

2 . 03223 

TLOC* 

1 .00000 

XMA* 

55.12334 

XM3* 

10.00000 

XMU* 

30.87615 

KR*  2 

i  * 

^  r*  ro"1/ or 

1  a  K/AjrJ  ^  O 

TOTT* 

0 . 00000 

TLOC  * 

2 . 22020 

XMA* 

50.21920 

XMS* 

10.00000 

XMW* 

39.73080 
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Example  shown  in  Figure  9: 
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shown 

in  Figure  10. 
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